Shock waves in disordered media 



N. Ghofraniha^^, S. Gentilini^, V. Folli^, E. DelRe^'^ and C. Conti^'^ 
^ IPCF-CNR c/o Dipartimento di Fisica - Universitd La Sapienza, P. A. Moro 2, 00185, Roma, Raly 
^ ISC- CNR c/o Dipartimento di Fisica - Universitd La Sapienza, P. A. Moro 2, 00185, Roma, Raly 
^ Dipartimento di Fisica - Universitd La Sapienza, P. A. Moro 2, 00185, Roma, Raly 
Corresponding author: neda.ghofraniha@romal.infn.it 
(Dated: November 13, 2012) 

We experimentally investigate the interplay between spatial shock waves and the degree of disorder 
during nonlinear optical propagation in a thermal defocusing medium. We characterize the way the 
shock point is affected by the amount of disorder and scales with wave amplitude. Evidence for 
the existence of a phase diagram in terms of nonlinearity and amount of randomness is reported. 
The results are in quantitative agreement with a theoretical approach based on the hydrodynamic 
approximation. 



PACS numbers: 

Laser beams propagating in nonlinear media undergo 
severe distortions as the power is increased: spreading 
due to diffraction can be progressively reduced through 
self narrowing, up to the generation of solitons [1, 2 and 
dissipative and dispersive shock waves (SWs) [3-8 , thus 
fostering the formation of a variety of nonlinear waves. 
The way these are affected by disorder is a leading main- 
stream of modern research [91412]. Attention is given to 
the competition between strongly nonlinear and coher- 
ent phenomena, and their frustration due to random- 
ness and scattering; recent theoretical investigations deal 
with general frameworks described by "phase-diagrams" 
in terms of the two parameters: the amount of nonlinear- 
ity and of disorder [13]. However, no direct experimental 
nonlinearity-disorder phase diagram has been reported. 

The case of SWs is specifically relevant [T4UI6] , as 
they represent a strongly nonlinear and coherent oscil- 
lation (the undular bore) [H [T7H2Q^ and are expected to 
be strongly affected (and eventually inhibited) by disor- 
der, at variance, e.g., with solitons, which can survive a 
certain amount of randomness (see, e.g., [211 [22]). This 
leads to the direct opposition between the two effects: 
on one hand increasing the nonlinearity favors the shock 
formation, on the other hand disorder-induced scattering 
limits this phenomenon. This is relevant in colloidal sys- 
tems [HI [23H27] where disorder is unavoidable, as well as 
in out-of-equilibrium photorefractive nonlinearities [28 , 
optical fibers [5l [29] , and also in Bose-Einstein condensa- 
tion [3QH32] and acoustics [33 . 

In this Letter, we report on the direct experimental ev- 
idence of the competition between SWs and disorder, and 
support our experiments by a theoretical model based on 
the hydrodynamical approximation. We measure the first 
phase-diagram (where the order parameter is the posi- 
tion of the formation of the shock) for nonlinear waves in 
terms of disorder and nonlinearity, and characterize the 
scaling laws for the random SWs formation and propa- 
gation. 

Experiment — We use dispersions of silica spheres of di- 
ameter 1/im in O.lmM aqueous solutions of RhodamineB 
displaying a thermal defocusing effect due to partial light 



absorption [TH \TE\ [34l [35] . To vary the degree of disorder 
several silica concentrations c are prepared, ranging from 
0.005 w/w to 0.03 w/w, in units of weight of silica parti- 
cles over suspension weight. A continuous- wave laser at 
wavelength A=532 nm is focused on the input facet of the 
sample (beam waist 10/im). To detect light transmit- 
ted at the exit of the samples, the aqueous solutions are 
put in 1mm X 1cm x 3cm glass cells with propagation along 
the 1mm vertical direction (parallel to gravity) to mod- 
erate the effect of heat convection. Transverse images of 
the beam intensity distributions are collected through an 
objective and recorded by a 1024x1392 pixel CCD cam- 
era. All measurements are performed after the temper- 
ature gradient reaches the stationary state and the par- 
ticle suspensions are completely homogeneous. The loss 
mechanisms in our samples are absorption and scatter- 
ing. The measured loss length (absorption + scattering) 
is L 1.6mm for pure dye solution and L 1.2mm for 
the sample with the highest concentration c (and hence 
highest losses) . These values are obtained by fitting with 
exponential decay the beam intensity versus propagation 
distance Z. The fact that L is always greater than the 
position of the shock point Zg (measured below) allows 
us to neglect losses at a first approximation in our theory. 
In addition, we find that the scattering mean free path is 
of the order of millimeters for all the considered samples. 
In Fig. [l] we show images of the transmitted beam (on the 
X-Y transverse plane) for different input laser powers P 
and various concentrations c. The profiles display post- 
shock rings with outer rings being more intense than the 
inner ones, as typical for dispersive SWs from Gaussian 
beams [14 . The number and the visibility of the oscil- 
lations increase with P and decrease with c, evidencing 
that SWs are sustained by nonlinearity and inhibited by 
disorder. 

We then investigate the onset of the shock along the beam 
propagation direction: we use 1cm x 1cm x 3cm glass cells 
(propagation along 1cm), and top images are collected 
through a microscope and recorded by the CCD cam- 
era. The effect of disorder on SWs along beam prop- 
agation is reported in Fig |2| where collected images of 
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FIG. 1: (Color online) Images of transmitted intensity for 
different input power P and particle concentration c: a) P=5 
mW, c= 0, b) P=400 mW, c= 0, c) P=5 mW, c=0.017 w/w, 
d) P=400 mW, c= 0.017 w/w e) P=5 mW, c= 0.030 w/w, 
f) P=400 mW, c=0.030 w/w. Superimposed curves show the 
measured section of the intensity profiles. 



does not have any maximum in the observation window 
Lo — 1 mm. In this regime Zs is not the position of the 
peak of the steepness but that of the highest steepness 
available in the observable Z range. We define the value 
of P at which Zs starts to decrease as the treshold power 
between shock and non-shock regimes and map the phase 
diagram in Fig. [3b. 
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FIG. 3: (Color online) (a) Measured shock point Zs Vs P for 
various c; (b) disorder-power phase diagram with shock and 
non-shock regimes obtained from the data in panel (a): dots 
correspond to the threshold powers, dashed line and the dot- 
dashed line are the boundaries as estimated by the theory. 



the transverse distribution of the beam intensity versus 
Z for different input power P and silica concentrations c 
are shown. In accord with Figjl] shock inhibition by dis- 
order in Fig [2] is evidenced by the reduction of the beam 
aperture and the disappearance of the undular bores. 





Disorder 








(a) 


(0 


(e) 











CL 








cu 

LO 








_l 


(b) 


(d) 


(f) 


y 


—150 ^im 







FIG. 2: (Color online) Beam propagation as observed from 
top fluorescence emission for different input power P and par- 
ticle concentration c: a) P=8 mW, c= 0, b) P=450 mW, c= 
0, c) P=8 mW, c=0.017 w/w, d) P=450 mW, c=0.017 w/w 
e) P=8 mW, c=0.030 w/w, f) P=450 mW, c=0.030 w/w. 

We identify the point of shock formation Zs as the Z 
corresponding to the maximum of steepness of the inten- 
sity profiles (details in [16 ). We follow this procedure 
for all images and we report in Fig. [3^ Zs versus P for 
different concentrations c. Two effects are evident: 1) 
for increasing power P, Zs decreases, corresponding to 
the speed-up of the shock formation caused by the aug- 
mented nonlinearity; 2) for increasing concentration c, Zs 
increases, as disorder delays shock formation up to its to- 
tal cancellation observed for c=0.030 w/w (star symbols). 
The plateau at low P indicates that shock is not occur- 
ring, as the steepness of the profiles increases with Z but 



Theory — In order to theoretically analyze the exper- 
imental results, following previously reported analyses 
O [22, we use the hydrodynamic approximation. We 
start from the paraxial wave equation for the field com- 
plex envelope A, in presence of disordered local Kerr 
medium with refractive index perturbation An = 77,2/ + 
AnR{X,Y,Z) to the bulk index no, with / = |Ap 
the optical intensity, n2 < the Kerr coefficient and 
Ani^(X, F, Z) a random perturbation, 

2ik^ + VlyA + 2k^—A = 0. (1) 
oZ ' no 

In Eq. ([1]), /c = 27rno/A, and we neglect spatial nonlo- 
cality and losses, because they do not qualitatively affect 
the scenario (as will be reported elsewhere). The corre- 
sponding dimensionless equation for the normalized field 
t/j = A/v^, with Iq the input peak intensity, is 

ie^^'^Vly^-\^\'i^^UR^ = 0, (2) 

where {x,y) = {X,Y)/wo, z Z/L, L y^L^L^i, 

e \/Lni/Ld, Ld kwl, Lni no/(/c|n2|/o), and 
Ur = AnR/n2lo is the ratio between the random in- 
dex perturbation and the nonhnear one {wq = 10/im). 
Because of the huge number of particles randomly dis- 
tributed within the optical beam and the low index con- 
trast (the refractive index is 1.46 for silica and 1.33 for 
water), nR{x^y^z) can be taken as a random dielectric 
noise mainly acting on the phase of the propagating 
beam; below we show that such an assumption allows 
to retrieve theoretical results in quantitative agreement 
with experiments. In the hydrodynamic limit e ^ 
{Lni << Ld)^ the propagation of the field intensity can 
be separated by that of the beam phase, this results into 
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the equation of motion for the phase chirp identical to 
that of a unitary mass particle (due to the cylindrical 
symmetry of the system we limit to the x-z variables): 

(fx dU dUn dU 
dz'^ dx dx dx 

In Eq. ([3| /7 = exp{—x'^/2) is the deterministic po- 
tential from the nonlinear part due the Gaussian beam 
profile. r]R = —dUn/dx is taken as a Langevin force 
that we assume with Gaussian distribution, such that 
{Vr{^)Vr{^')) = r]'^S{z — z')^ with strength of disorder 
measured by 7/ = {{dUn/dr)'^)'^^'^ =< An^ >^/^ /|^2|^o, 
and the brackets denoting statistical average. We stress 
that in the following we solve eq.(|3| for several values 
of X taking for each of them an independent realization 
of the noise r]R{z)^ this allows to neglect the dependence 
of rjR on X. We stress that for reasons of symmetry the 
disorder in the two transverse directions are independent. 
The simplest and effective theoretical approach is to con- 
sider a one-dimensional reduction. Because of the disor- 
der averaging cylindrical symmetry is preserved, as also 
experimentally demonstrated. 

In Figj4^,b, we show several of these trajectories result- 
ing from initial uniformly distributed position in the x 
axis and zero initial velocity v = dx/dz, as obtained by a 
stochastic Runge-Kutta algorithm [36]. Upon propaga- 
tion the particles collide and, in the absence of disorder, 
the shock is signaled by the coalescence of multiple tra- 
jectories (Fig. [4^); in the phase space of v and x (FigJIJ^), 
these correspond to the folding of the velocity profile into 
a multivalued function when increasing z, which induces 
the wave-breaking phenomenon. In the presence of dis- 
order, the particles tend to diffuse, as evident from their 
trajectories (Fig|4]3) and in the phase space (Fig|4]i); cor- 
respondingly, the propagation distance before their col- 
lisions is greater for their random walk and the shock is 
delayed in the direction. 

Figures |5^,b show the numerically obtained histograms 
of the particle positions at various propagation distances. 
If compared with the ordered case in FigjS^, disorder in- 
duces a spreading of the particle distribution. We extract 
the position of the shock Zs = Zs{r]) as that approxi- 
mately corresponding to the maximum of the histogram 
(precisely, as the mean value among the positions for 
which the histogram is above the 90% of its maximum, 
to limit fluctuations). This allows us to determine Zg 
for various amounts of disorder r] (in the ordered case 
Zs{0) = 2.5). Figure [s]^ shows Zs{r]) vs. disorder degree 
for 10^ particles and demonstrates that the shock process 
is delayed by disorder. 

As discussed above, in the absence of disorder, shock 
appears in the experiments only above a threshold power 
(see FigjS^i): from a theoretical point of view this thresh- 
old arises from the fact that the hydrodynamic model 
{Lni « corresponding to e ^ 0) is valid only at 
high nonlinearity, hence no shock is expected at low 
power. Moreover, in the hydrodynamic limit the po- 
sition z = Z ^/\rl2\PJ7T / Wq and the shock position Zg 



scales as as experimentally investigated in [16]. 

Following our theoretical approach, this scaling is 
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FIG. 4: (Color online) Trajectories of colliding particles form- 
ing shock versus z: (a) without disorder (77 = 0) and, (b), for 
?7=0.1; (c) shock profile in the phase space for 77 = and, (d), 
for ?7=0.1 {z varies in the range [0,3]). 
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FIG. 5: (Color online) Theoretical histograms of particle po- 
sitions for ?7 = (ordered case, panel a) and for 77 = 0.2 
(panel b) ; (c) theoretical normalized shock position Zs versus 
amount of disorder 77 (black continuous line) and comparison 
with the measured Zs Vs concentration c (red squares). 

maintained in the disordered case and the shock posi- 
tion is delayed, such that Zs{r]) = Zs{r])wQ/ ^\n2\P/T^^ 
Zs is compared in FigjS]^ with our experimental results 
(|n2| = 2 X 10~^^m^W~^) revealing quantitative agree- 
ment; discrepancies between experimental and theoret- 
ical Zs are ascribed to the several adopted simplifying 
assumptions in the latter. 

In addition, from the theory another threshold arises, 
corresponding to the existence of a critical value for 
the amount of randomness above which no shock is ob- 
tained. Indeed disorder becomes dominant with respect 
to nonlinearity when Ur \^ greater than the determinis- 
tic part [/, such that the hydrodynamic model in Eq.(|3| 
is no longer valid. More precisely, Eq.([3| holds true 
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as long as 7^ ^ 1, above this value no shock is ex- 
pected. As U is of the order of unity, this corresponds 
to = {Aniy/y\n2\Io ^ 1, with {Aniy/^ = 

cpH2o{nsi02 - ^H^o)/ PSi02^ being nsi02 {nH^o) and 
Psi02 {PH20) the refractive index and the density of the 
Si02 {H20)^ respectively. That is, no shock is expected 
when the random index perturbation Anji becomes com- 
parable with the nonlinear perturbation 712/, such that 
material fluctuations are so pronounced that the non- 
linear effect is totally masked. In our experiments we 
have |n2|/o — 10~^, and we have that the condition 
(An^)^/^ = |n2|/o — 10~^ is found for a concentra- 
tion of c = 0.03 w/w (see figure [sj^). In addition, the 
other boundary line in the phase diagram scales as a/P 
(dot-dashed in figure [s]^), as retrieved from the condition 

Conclusions — By the direct visualization of a laser beam 
propagating in a liquid random nonlinear system and by 
imaging the transmitted light at the exit of the samples, 
we show that the gradual increase of disorder hampers 



shock wave formation, up to its total inhibition. Such 
a transition has been quantitatively characterized and 
results into the first direct measurement of the phase- 
diagram of nonlinear propagation in terms of disorder 
and nonlinearity, as also supported by a theoretical model 
based on the hydrodynamic approach. These experi- 
ments open the way to further investigations concerning 
the interplay between disorder and nonlinearity, as the 
identification of glassy and superfluid phases of light and 
related phenomena. 
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